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AN ANALYSIS ОЕ LAMINAR FREE-CONVECTION FLOW AND HEAT TRANSFER 
ABOUT A FLAT PLATE PARALLEL TO THE DIRECTION 
OF THE GENERATING BODY FORCE 1 


By ЁЯмом OsrRACH 


SUMMARY 


The free-convection flow and. heat transfer (generated. by თ 
body force) about a flat plate parallel to the direction of the body 
force are formally analyzed. and. the type of flow 18 found. to be 
dependent on the Grashof number alone. For large Grashof 
numbers (which are of interest in aeronautics), the flow 18 of 
the boundary-layer type and the problem is reduced in а formal 
manner, which is analogous to Prandtl’s forced-flow boundary- 
layer theory, to the simultaneous solution of two ordinary 
differential equations subject to the proper boundary conditions. 

Velocity and temperature distributions for Prandtl numbers 
of 0.01, 0.72, 0.783, 1, 2, 10, 100, and 1000 are computed, and 
it is shown that velocities and Nusselt numbers of the order of 
magnitude of those encountered in forced-convection flows may 
be obtained in free-convection flows. The theoretical and 
experimental velocity and temperature distributions are in good 
agreement. 

A flow and თ heat-transfer parameter, from which the impor- 
tant physical quantities such as shear stress and heat-transfer 
rate can be computed, are derived as functions of Prandtl num- 
ber alone. Comparison of theoretically computed values of the 
heat-transfer parameter with values obtained from an approzi- 
mate calculation and experiments yielded good agreement over a 
large range of Prandtl number. Agreement between the 
theoretical values and those obtained from a frequently used 
semiempirical heat-transfer law was good only in restricted 
Prandtl number ranges (depending on an arbitrary constant). 


INTRODUCTION = 


Two important types of fluid flow problems involving heat 
transfer are those of forced and those of free convection. 
By forced-convection flow is meant flows maintained mechan- 
ically as, for example, by a pressure drop or an agitator. 
Tree-convection flow, on the other hand, results from the 
action of body forces on the fluid, that is, forces which are 
proportional to the mass or the density of the fluid. The 
flow is generally produced in the following manner: Consider, 
for example, a fixed object (such as a plate) in a quiescent 
fluid subject to a body force. When the plate is at the same 
temperature as the surrounding fluid, the body forces acting 
on the fluid are in equilibrium with the hydrostatic pressure 
and no flow ensues in the steady state. If a temperature 
gradient normal to the body force is imposed by heating 
(or cooling) the plate, there will exist a defect (or excess) of 


body force because of the decreased (or increased) density, 
with the fluid closer to the plate having the greater defect 
(or excess) than that away from the plate. ‘This unbalance 
of the forces causes the fluid to be accelerated with the 
particles nearer the plate moving more rapidly than those 
farther from the plate. Free-convection flow has usually 
been considered to be generated in a gravitational field where 
the previously mentioned defect or excess of body force was 
the Archemedian (buoyancy) force. However, since cen- 
trifugal forces are also proportional to the fluid density, free- 
convection flows can also be set up by the action of such 
forces. (See ref. 1 for a more explicit discussion of the 
development of. free-convection- flows by centrifugal forces.) 

Free-convection flows produced by centrifugal forces are 
now of practical importance in aeronautics because many 
aircraft propulsion systems contain components (such as gas 
turbines and helicopter ram jets) which rotate at high 
speeds and in which heat is being transferred. The method 
of free-convection cooling of gas-turbine rotor blades where 
the centrifugal forces create a free-convection flow of the 
coolant in the blade passages is an example of a practical 
application of the free-convection phenomenon in aero- 
nautics. Also, free-convection flow due to centrifugal force 
is superimposed on the flow through helicopter ram jets and 
on the flow of cooling air in hollow rotor blades of air-cooled 
turbines. and, under proper conditions, can appreciably 
influence the resultant flow and heat transfer. | 

‘As a simplification of the many free-convection problems 
which are now of some consequence in aeronautics, consid- 
eration is here given to the special case of free-convection 
flow about a flat plate parallel to the direction of the gener- 
ating body force. The experimental and theoretical consid- 
erations of Schmidt and Beckmann (ref. 2) concerning the · 
free-convection flow of air subject to the gravitational force 
about a vertical flat plate constitute the most complete 
treatment of this subject up to the present time. Eckert 
(ref. 1) as well as others has further verified and extended 
the experimental results of Schmidt and Beckmann, and 
Schuh (ref. 3) has extended the numerical calculations by 
computing the velocity and temperature distributions for 
several Prandtl numbers different from that for air. How- 
ever, all the theoretical work in these references is based on 
the incompressible equations in which the density (or tem- 
perature) variation is introduced in the buoyancy term alone. 


1 Supersedes NACA ТМ 2635, “Ап Analysis of Laminar Free-Convection Flow and Heat Transfer about a Flat Plate Parallel to the Direction of the Generating Body Force" by Simon 


Ostrach, 1952. 
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Various terms are omitted from the equations at the start 
on the basis of either intuitive arguments or no arguments 
at all. Although a theoretical development made in such а 
manner led to good final results, the significance of all the 
important factors associated with the free-convection flow 
phenomenon is not obtained from such an analysis. 

The problem of free-convection flow as produced by a body 
force about a flat plate in the direction of the body force was 
studied at the NACA Lewis laboratory during 1951 and is 
treated in a formal and more general manner herein. The 
method used is somewhat similar to that used in reference 4 
wherein consideration was given to the free-convection flow 
at high Grashof numbers in a horizontal cylinder which had a 
variable surface-temperature distribution. The application 
of this method to the present problem leads to a development 
which is analogous to Prandtl’s treatment of high Reynolds 
number forced-convection flows. Although the final equa- 
tions obtained by this method are the same as those of 
Schmidt and Beckmann, this more general approach not only 
clearly demonstrates the significance of all the important 
parameters and assumptions and hence leads to a better 
understanding of this type of flow but also indicates the 
quantitative limitations of the theory. In addition, the 
numerical solutions of references 2 and 3 are herein extended 
to cover & more complete range of parameters. The new 
calculations yield information on the free-convection flow for 
Prandtl numbers corresponding to those of liquid metals, 
gases, liquids, and very viscous fluids. 


ANALYSIS 
STATEMENT OF PROBLEM AND BASIC EQUATIONS 


The steady-state equations expressing the conservation of 
mass, momentum, and energy for a compressible, viscous, 
and heat-conducting fluid subject to a body force together 
with an equation of state govern the flow and associated 
temperature distribution about the plate. These equations 
in Cartesian tensor notation are (see ref. 5), respectively, 
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(А complete list of the symbols used herein is given in 
appendix A.) For the two-dimensional case, equations (1) 
to (4) represent а system of five equations in the five depend- 
ent variables Ui, Uz, p, P, and T. For later use, оно (4) 
can be written 

dp=p(K dP—B dT) (4a) 


where К and В are the coefficients of isothermal compressi- 
bility and volumetric expansion, respectively (see ref. 6). In 
addition to a general state equation, such as is given in 
equations (4) or (48), it will be convenient at times in the 
discussion to refer to some specific state equation. To this 
end, the'equation of state for an ideal gas 


P=pRkT (4b) 
will be used. 

Particular consideration is here given to the two-dimen- 
sional free-convection flow about а semi-infinite vertical flat 
plate. The X;-axis of the coordinate system is taken along 
the plate and the X;-axis, normal to it. No distinction is 
made as to the specific type of body force acting, for example, 
gravitational or centrifugal, but һе force is assumed to be 
acting in the vertical direction only (that is, parallel to the 
plate). Centrifugal and Coriolis forces which are connected 
with flows on curved paths and with rotating systems 
generally vary with position and velocity. However, in 
order not to make the analysis unduly complicated, the body 
force is taken to be constant. | 

In order to define the problem clearly, а choice must still 
be made of the position of the origin of the coordinate sys- 
tem. Before making a definite decision on this point, note 
that for constant plate temperatures there are four permu- 
tations of the body-force direction (either upward or down- 
ward) and the plate thermal condition (either heated or 
cooled) which will lead to free-convection flows. Once the 
position of the edge of the plate, which is also to be the 
origin of the coordinate system, is decided, there are two 
combinations of the body-force direction and plate thermal 
condition that will yield flows which proceed away from the 
edge. It is this type of flow that is amenable to the type of 
analysis to be made here. This point will be more fully dis- 
cussed subsequently. И the edge of the plate (recall that ე, 
semi-infinite plate has but one edge) is taken at the bottom of 
the plate (that is, the plate extends to + in the X,-direc- 
tion), the two combinations leading to flows in the proper 
direction (upward in this case) are, respectively, the body 
force acting downward with a heated plate and the body 
force acting upward with a cooled plate. The equations 
developed for one of the cases reduce directly to those for the 
other. The remaining two permutations, namely, the body 
force acting downward with a cooled plate and the body force 
acting upward with a heated plate, would yield flows which 
proceed downward or toward the edge of the plate if this 
edge were taken at the bottom of the plate. This type of 
flow would violate a physical condition of the problem which 
states that the flow starts at the plate edge. The latter 
combinations hence will not be considered further. 

Because the two acceptable configurations can be reduced 
essentially to one, for the development to be given here, the 
origin of the coordinate system will be taken at the bottom 
of a heated plate, with the body force acting downward. 
The assumption is now made that the viscosity and thermal- 
conductivity coefficients are functions of the Ca a 
alone and obey the following laws: 
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The choice of the body-force direction together with equa- 
tions (5) alters equations (2) and (3) so that they become 


(5) 
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Note that the only nonzero component of the body force is 
the JX;-component. 


BOUNDARY CONDITIONS 


The boundary conditions associated with the given prob- 
lem are that: 

(a) The fluid must adhere to the "A (the no-slip condi- 
tion of viscous flows) and the plate must be а streamline, 
or mathematically, 


— 


0. (<. 0)=0. (5,0) =0 | (8) 


(b) The temperature of the fluid at the plate must be equal 
to the plate temperature, that is, 


T(X,,0)— T (9) 


(c) The velocity U, at large distances from the plate must 
be undisturbed, or 


U(X," )=0 (10) 


(d) The temperature at large distances from the plate 
must be equal to the undisturbed fluid temperature, or 


T(X,9)— T. (11) 
SIMPLIFICATION OF EQUATIONS 
Let a small quantity e now be defined as 


e—B(To— To) (12) 


which is а measure of the magnitude of temperature variation 
in the flow field. The coefficient of volumetric expansion В 
is generally of the order of magnitude between 107? and 10+ 
(see table 15 of ref. 7, for example) and for gases, 8— 1/7. 
(Thus, for gases, if В is taken to be constant, «= (To— T«)/T« ; 
that is, e is the relative temperature difference.) The 
coefficient В will be assumed constant. Because in the 
Bteady state flow ensues only when there is à temperature 


821690—55———6 


-variation in the fluid, the free-convection velocity should 


then depend directly on e, and the variations in pressure and 
density (from the static, -=0, case) due to the temperature 
differences should also depend on e. Thus 


pss (551 =") us (13) 
P=P,+P eo (14) 
p=pstpn€y (15) 
T=T (1+ еб) (16) 


where —fx denotes the X,-component of the body force рег 
unit mass, и, c, თ, and 0 denote dimensionless functions 
(which, in general, can be functions of e), [13 some character- 
istic length (for example, the distance from the edge of the 
plate to the point of interest), P, and p, are the pressure and 
the density, respectively, for the static case (U;=0 or 
є=0), and P. and po denote constant values of the pressure 
and the density (that is, the values if no force field were 
present) defined by the state equation (in the case of a gas, 
in particular, Р.=р»ВТо.): Because there is no character- 
istic velocity associated with the type of flow under con- 
sideration, the velocity is dimensionalized by the factor 
given in parentheses on the right side of equation (13). 

In order to determine the static quantities, it will at first 
be convenient to consider the particular case of & gas. "The 
problem 13 then considered with the temperature uniform 
throughout the flow field &t the value Т. (therefore there. 
will be no flow and U;=0). For this situation, equations 
(4b) and (6) become 


P,—p,RT, (17) 
and | 
ont рг}х= 0 
(18) 
ОР, | 0 
O Xs | 


(It should be noted that eq. (18) expresses the physical fact 
previously stated that the body force and hydrostatic 
pressure are in equilibrium for the static case.) Substitution 
of equation (17) mto equation (18) leads to 


3 


and equation (19) together with equation (17) and the 
equation defining Р. and р. yields 


P,=P, exp (- ИЕ (19) 


Ps= p, OXP (2 X)-». exp (-5e- X) (20) 


If the exponential in equation (20) is expressed in terms of 
its series expansion, that equation becomes 


P= р. (1-25 Xu. ) (21) 


66 REPORT 1111—NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


À computation of the second term in the parentheses of 
equation (21) for the case of air under normal conditions 
with fy=g and the fact that X, is of the order of magnitude 
| show that р-ё ИР. -10-*1/foot. For the type of problem 
under consideration, } will always be of unit order of magni- 
tude so that even if the body force fx represents a cen- 
trifugal force many times that of gravity, the inequality 
р-јх Х1Р..<1 may still be satisfied. Thus, in the subsequent 
development it will be assumed that p,cpo.. ‘This assump- 
tion, which was justified by the computation for the case of 
& gas, is expected to be reasonable for other fluids as well. 
'The physical interpretation of this assumption is that under 
static conditions (є==0), the density (or pressure) is not 
affected by the force field. 

In order that all quantities in the following equations be 
dimensionless, it is further necessary to define z,— Хі, 
where the z, are now dimensionless space coordinates. Sub- 
stituting these new coordinates along with equations (13) to 
. (16) into equations (1), (6), (7), and (4а) and noting equations 
(18) and that p,cp. yield, on neglection of terms of higher 
order in e compared with those of order e, 
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where P../p.fxl=NGr and the Grashof number Gr and the 
Prandtl number Pr are defined as 


9 f. 
Po Fx 
a 


and 


Рта = 


ae 


Physically speaking, the Grashof number represents the 
ratio of the body forces to the viscous forces. The free- 
convection flows of interest here are those associated with 
large Grashof numbers. The factors К and В in equation 
(26) may well be taken to be constants (see ref. 6). 

The boundary conditions (eqs. (8) to (11)) in terms of the 
new dimensionless variables are 


и (21,0) —5(2,,0) =0 (27) 


Об) — 575 (28) 


пл (21, ) —0 (29) 


0(21,9% ) =0 (30) 
Thus, to a first approximation, equations (22) to (26) to- 
gether with the boundary conditions replace the original 
equations and boundary conditions. (Note that for gases, 
BI.—l) 

The prime assumption made in this analysis is that the 
higher-order terms in e are negligible, which implies that є 
is small, and consequently, that the temperature difference 
ог В is moderately small. It is a consequence of this assump- 
tion alone that the basic equations were simplified to equa- 
tions (22) to (26). wherein the viscosity term in the energy 
equation is neglected and the only coupling of the momentum 
and energy equations occurs by means of the body-force | 
term in equation (23). As a result of this assumption, the 
variations of the viscosity and heat-conductivity coefficients 
with temperature are also negligible. Without any discus- 
sion, the authors of reference 2 start directly from simplified 
equations of the same form wherein the pressure terms in the 
energy equation were also neglected. In reference 3 some 
intuitive arguments are given to justify the simplified 
equations. 

It is now convenient to revert to the more familiar notation 
where 2=21, 2/= თეა, U=U, and =u. Equation (22) implies 
the existence of a stream function V such that 


and 


(31) 


where subscripts denote differentiation. Applying equation 
(31) to equations (23) to (25) yields, respectively, 


az Ah A= иф vel toe (32) 
% (Ay) EM уу: Viu TN Oy (3 3) 
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The boundary conditions (eqs. (27) to (30)) become 
Vs (2,0) = №, (2,0)— 0 (85) 
1 
(x, 07—57 (36) 
(37) 


p(z, œ)=0 (z, о )= 0) 


Equations (32) to (34) and equation (26) form the system 
of equations for the four unknown functions V, 0, o, and с 
of the problem. The system is nonlinear, and therefore 
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further simplification of the equations would be desirable. 
Just as in the case of forced-convection flows where the 
Reynolds number determines the type of flow or, in mathe- 
matical terms, the type of solution, the Grashof number is 
the prime factor for free-convection flows. For the case of 
small Grashof number, it can be seen from equations (32) 
to (34) that а perturbation in the small parameter Gr will 
yield а, system of linear equations. For Grashof numbers 
of unit order of magnitude, no further important simplifica- 
tion сам be made and the solutions would have to be obtained 
numerically. For the other limiting case, that of large 
Grashof numbers (which 15. the case under consideration 
herein), it would, at first thought, appear that some simpli- 
fication could be obtained by performing & perturbation in 
the small parameter 1/Gr. However, this would then 
imply that the term containing the highest-order derivatives 
(the left term in equations (32) to (34)) could, among others, 
be neglected. (This argument would also imply that the 
body-force term ф in equation (32), which is essentially, 
causing the flow, could also be neglected.) The omission 
of the highest-order derivatives from consideration, however, 
would lead to solutions which would not satisfy all the 
boundary conditions. Problems of this type are referred 
to as singular perturbation problems. For further dis- 
’ cussions of singular perturbation problems, see references 
8 and 9. 

Equations in which a small parameter multiplies the 
highest-order terms are said to be of the boundary-layer 
type, because in order for solutions which satisfy all the 
boundary conditions to be obtained, the highest-order terms 
must be considered near the boundary. This fact implies 
the existence of a thin region, called the boundary layer, 
wherein the functions vary rapidly from the value at the 
boundary to that in the flow outside this layer. The con- 
clusion to be drawn from the preceding discussion is that 
for large Grashof numbers the flow is of the boundary- 
layer type. Schmidt and Beckmann (ref. 2) also made 
the boundary-layer assumptions in their theoretical develop- 
ment, and these assumptions were justified on the basis of 
their experimental observations. The Grashof numbers for 
their experiments were of the order of 8Х 10°. 

In view of the fact, previously discussed, that highest- 
order derivatives of each dependent variable as well as of 
those terms of physical importance (as, for example, the 
body-force term) must be retained in the boundary 1 ს it 
is convenient to make both sides of each of the equations 
of the same order in Gr. In this way, as will be shown, the 
equations will be further simplified. It is thus convenient 
to make the following transformations in the system of 
equations (32) to (84) and (26) and then to retain only the 
dominant parts (that is, those multiplied by Gr to the 
highest power) of each individual term. 

Let j—Gr'y, y-Gr*y, o—Gr'v, ф=ф, and 0—0. 
equations (32) to (34) and (26) become 


Then 
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[6 now can be seen that by proper choice of г, 8, and £ а 
transformation of the type given provides & means for 
making the important terms in the differential equations of 
the same order in Gr. Thus if r—X, s=—%, and t=—1, 
equations (38) to (41) become 


Vos 9— Ves — Vd + No. (42) 
AX. 

No,—0 (ar 4) eo (43) 

60;=Рг(4.0,— 9.8; (44) 

45--ВТ.46=0 (45) 


More generally, if N is very much different from unit order 
of magnitude, а value of ¢ can always be chosen (depending 
on AN) such that equations (42) to (45) are obtained. (For 
any negative ¢ less than —1, the last term of eq. (42) 
will also disappear.) | 

There are now several important points to be discussed 
concerning the transformation just made and the resulting 
simplified equations. First, it should be noted that the 
transformation is merely & formal expression of the bound- 
ary-layer assumptions first made by Prandtl and hence the 
solutions will be asymptotic for large Gr. Second, the 
second equation of motion here also reduces to state that 
the pressure across the boundary layer is constant. Third, 
the pressure terms in the energy and state equations are 
here found to be negligible. This fact verifies a priori 
assumptions made by others from the physics of the problem. 
Finally, note that integration of the general state equation 
(independent of pressure) as now given by equation (45) 
leads to 

| oT 8T.6—0 (46) 
where the constant of integration has been taken as zero 
without any loss of generality. For the particular case of 
a gas, B=1/T. so that equation (46) becomes 


o+6=0 


The boundary conditions (eqs. (35) to (37)) now can be 
written 


‚ (2,0)=¥,(2,0)=0 (47) 
- 1 
9 (2,0) gm- (48) . 
Vs (2, ©)=8(2, )—0 (49) 


If now it is assumed that 7,—0 in equation (42) since con- 
sideration is here being given to a flat plate, and if დ is 
eliminated from equation (42) by use of equation (46), 
there results the system of equations 
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Yn -ВТ.8= 9—4 = (50) 
05==Рт (459, — 1,05) (51) 


Thus the problem has been reduced to the solution of the 
two simultaneous partial differential equations (eqs. (50) and 
(51)) subject to the boundary conditions (eqs. (47) to (49)). 

Final simplification of the equations is made by applica- 
tion of the so-called similarity transformation of boundary- 
layer theory. Thus, let 


=>” (59) 
(4а) 
and 
V—(42)* F(x) (53) 
z На) 
р (54) 


Then equations (50) and (51) are reduced to the following 
ordinary differential equations: 


F”! 43 FF” —2F"°+H=0 (55) 


H'' ОР F H'=0 | (56) 


where the primes denote differentiations with respect to 7. 
The reciprocal one-fourth power similarity as given in equa- 
tion (52) is characteristic of free-convection flows just as the 
reciprocal square-root type is characteristic of the forced- 
convection flows. "The boundary conditions become 


JF" (0) =Р(0) =0 (57) 
H(0)—1 (58) 
F'(*)—H(9)—0 (59) 


The use of a transformation like equation (52) essentially 
specifies an additional boundary condition, namely, that the 
conditions to be satisfied at у= თ (or n= œ) should also be 
satisfied at х=0. It is for this reason that the flows pre- 
viously discussed which would flow toward the edge (down- 
ward) are not amenable to this type of analysis, for such flows 
would violate this additional condition, which essentially 
states that the boundary-layer development starts at the 
. edge of the plate. 


SOLUTION OF THE BOUNDARY-VALUE PROBLEM 


The solutions of the simplified equations (55) and (56), ' 


satisfying the boundary conditions as given by equations 
(57) to (59), were obtained by use of ап IBM. Card-Pro- 
grammed Electronic Calculator. A detailed account of the 
procedure followed in the determination of the unknown 
functions is presented in appendix B by Dr. Lynn U. Albers. 
The functions F and H together with their derivatives are 
given in table I for Prandtl numbers of 0.01, 0.72, 0.733, 1, 2, 
10, 100, and 1000. Even though the Prandtl number for air 
is taken as 0.72 in this report, the solutions for Pr—0.733 
were 9150 computed and are presented as а check with the 


Schmidt-Beckmann calculations wherein the value of Prandtl 
number of 0.733 was used. The particular values of the 
Prandtl numbers given were chosen to correspond to those 
for liquid metals, gases, liquids (such as water and oil), and 
very viscous liquids (such as glycerin or oils at very low 
temperatures). | | 
RESULTS 
VELOCITY AND TEMPERATURE DISTRIBUTIONS 


By means of the various transformations made in the 
analysis it can easily be verified that 


Us 
U 
o laf mex =” 
2./В(Т,—Т. Хо FOTE (я) (60) 
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(61) 
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Equations (60) to (62) relate the physical quantities to the 
dimensionless functions Ё and Н which are now known. 
The dimensionless velocity and temperature distributions as 
given by equations (60) and (61) are presented in figures 1 
and 2, respectively, as functions of 7 for the various values 
of Prandtl number. The computations made here agree 
with those for Pr—0.733 as given in reference 2 up to the 
third significant figure. For Pr—10, 100, and 1000, the 
present results agree in general with those of reference 3. 
Since only curves are presented in reference 3, the precision. 
of the agreement cannot be stated. 

The maximum values of the dimensionless velocity dis- 
tributions occur at larger values of the argument უ as the 
Prandtl number decreases and the velocities decrease with 
increasing Pr. It should also be noted that the dynamic 
and thermal boundary-layer thicknesses can be estimated 
from the abscissas of figures 1 and 2, respectively, and that 
for Pr>1 the velocity boundary layer is much thicker than 
the thermal boundary layer. | 

The occurrence of fx (ог Grx as given by eq. (60)), which 
may be very large for flows generated by centrifugal forces, 
in the denominator of the ordinate implies that velocities of 
appreciable magnitude can be associated with such free- 
convection flows. In particular, if ]х==108 feet per second 
squared, which is a reasonable conservative figure for pres- 
ent-day rotating systems, e—0.2 (which is within the limits 
cf the theory presented herein), and arbitrarily А=0.25 
foot, then the maximum velocity attained at a Prandtl 
number of 0.72 is approximately 125 feet per second. Tbis 
value of the maximum velocity could, of course, be doubled 
or even tripled under the proper conditions. One limitation 
to a calculation of this sort, аз can be seen by comparison of 
the denominators of the left and middle terms of equation 
(60), should be kept in mind; namely, the limiting Grashof 
number for laminar flows. In lieu of a complete stability 
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analysis on this type of flow, this limiting value is taken to be 
10°, as indicated in reference 10. Consideration of this lim- 
itation then .implies (see eq. (60)) that for large laminar 
velocities either va must be large or X must be small. 


COMPARISON WITH EXPERIMENTS 


Careful experiments of free-convection flows (as generated 
by gravitational forces) about vertical flat plates were made 
by Schmidt and Beckmann (ref. 2) in which velocity measure- 
ments at various points along the plate were made by means 
of a quartz-filament anemometer and the temperature meas- 
urements were obtained by means of manganese-constantan 
thermocouples. Eckert (ref. 1) performed similar experi- 
ments in which the measurements were made by means of & 


FIGURE ].—Dimenstonless velocity distributions for various Prandtl numbers. 


Zehnder-Mach interferometer. The results of both sets of ` 
experiments are in 2004 agreement, but since the data pre- 
sented in reference 2 by Schmidt and Beckmann appear in 
more detail, these data will be used for comparison with the 
theory. - 

The experiments of reference 2 were performed on two 
different (in that the edges were smoothed either symmetri- 
cally or not) 12- by 25-centimeter plates and on one 50- by 
50-centimeter plate. It should here be pointed out that the 
results for the two smaller plates were almost identical and 
that the flow was entirely laminar except near the outer 
edge of the boundary layer where the slight turbulence of 
the room air disturbed the measurements somewhat. (This 
effect was also observed by Eckert.) Large periodic oscilla- 
tions of the flow near the downstream edge of the larger plate 
were observed in addition to the slight turbulence near the — 
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FiauRE 2.—Dimensionless temperature distributions for various Prandt! numbers. 
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outer edge of the boundary layer. Hence the data from the 
larger plate should not be expected to yield completely satis- 
factory agreement with the laminar theory as presented 
here. 

Since the physical quantities can be expressed in terms of 
9, single variable as in equations (60) and (61), it 13 to be ex- 
pected that the data taken at the various points along the 
plates should all lie on a single line if the data are correlated 
according to equations (60) and (61). Thus for the smaller 
plates where (7-7.)=95.22° R and Т.=518.68° В, equa- 
tions (60) to (62) become  . 


=F" (n) 


U 
———— 63 
4.8624 X (68) 
T—518.68 
—095.99 ^" =H(n) (64) 
Y 
–. (65) 


The velocity and temperature distributions are so plotted in 
figures 3 ара 4, respectively, as are the curves computed 
theoretically for Pr=0.72. It can Бе seen that the agreement 
is in general very good for small values of 7 and somewhat 
less satisfactory though still rather good for the larger values 
of y. The scatter in the range of the larger values of 7 is 
believed to be caused by the previously discussed room 
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FIGURE 3.—Oomparison of small plate experimental and theoretical velocity distributions 
for Prandtl number of 0.72. 
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Етасве 4.—Comparison of small plate experimental and uo eee temperature distribu- 
tions for Prandtl] number of 0 


turbulence. It should also be noted that the points farthest 
away from the theoretical are those measured near tho 
leading edge. These points should not, of course, be expected 
to agree too well with the theory since the boundary-layer 


'| assumptions made in the theoretical development imply that 


the distance along the plate is large as compared with the 
‘boundary-layer thickness. Hence, this assumption is invalid 
near the leading edge. Schmidt and Beckmann obtained 
closer agreement between the theory and the experiments for 
the temperature data and poorer agreement for the velocity 
data by basing the kinematic viscosity coefficient in equation 
(62) on the plate temperature rather than on the undisturbed 


.| stream temperature as was done here. 


For the larger plate, (To— Т.) =83.7° К and To — 6527.14? 
R, so that equations (60) to (62) become 


U 
Ор 66 
L5224*X ^ (80) 
T— 527.14 
837 Hl) (67) 
1—83.93 >, (68) 
= | 


The velocity and temperature distributions for this experi- 
ment are plotted in figures 5 and 6, respectively, and again 
the theoretical curves for Pr=0.72 are included. In figure 5 it 
can be seen that for large 7 the agreement is rather poor, 
particularly for the data for both small and large values of X. 
The poor agreement for small values of X is again due to the 
theory limitation near the edge of the plate and for large 
values of X, to the fact that the flow was becoming turbulent 
"there. 
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On the basis of а simplified theory (that is, Бу use of 


tegrated momentum and energy equations and assumed 


FREE-CONVECTION FLOW AND HEAT TRANSFER ABOUT A FLAT PLATE 


4. 


0.718(Pr)* 
(0.9524-Pr)* 
(83.93) У/Х 


1) 
FIGURE 6.—Comparison of large plate experimental and theoretical temperature distribu- 


Чопа for Prandtl number of 0.72. 
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FIGuRE 7.—Dimexsionless flow parameter as function of Prandtl number. 
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FIGURE 8.—Dimenslonless heat-transfer parameter as function of Prandtl number. 
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velocity and temperature distributions), Eckert (see p. 162 


of ref. 7) obtained the approximate relation 


in 


for Prandtl number of 0.72. 
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The curve representing this equation is also presented in 
figure 8, and it closely approximates (to within about 10 
percent) the curve determined by the more exact considera- 
tions of this report over the entire Prandtl number range. A 
semiempirical equation as given in reference 11 relating the 
average (over the length X) Nusselt number to the Prandtl 
and Grashof numbers which has been used in the heat- 
transfer calculations up to the present is 


Nu,,—0.548 [(Pr)@r]* 


The constant 0.548 pertains specifically to air; for oil it 
should be 0.555 (see ref. 12) and for mercury, approximately 
0.33 (see ref. 13). In order to obtain local values from the 
average ones given by the last equation, it is merely necessary 
to multiply the average values by 0.75. (The determination 
of this reduction factor of 0.75 is discussed in appendix D.) 
Thus in terms of the local quantities the semiempirical 
relation becomes 


Nu=0.411 [(Р^) ах 
or 
Nu 


— = 0.581 (Prt 
(@тх/4) 
The curve given by this equation is also presented in figure 
8 and the agreement with the theoretical curve is very good 
for Prandtl numbers near unity, not so good for large Prandtl 
numbers, and very poor for the small Prandtl numbers. Of 
course, changes of the constants in the semiempirical relation 
ав previously discussed for the large or small Prandtl number 
cases (oil and mercury, respectively, for example) would 
cause the semiempirical curve to approximate the theoretical 
curve more closely. "The values of the heat-transfer param- 
eter obtained experimentally for mercury (Pr=0.03), air 
(Pr=0.72), water (Pr=7), and oil (Pr— 75.5, 115, 190, 224, 
275, 318, 308, and 442) are here reduced by the factor 0.75 
from the average values reported. The value for mercury 
is an average taken of four readings from a curve, since this 
experiment was the only one not reported in tabular form. 
From figure 8 it can be seen that all of the experimental 
values except those for the oil experiments are 10 very good 
agreement with the theoretically computed values. The 


data from the oil experiments, though not so good, show 
reasonable agreement (maximum error of approximately 20 
percent) with the theoretical curve and good agreement, as 
is to be expected, with the semiempirical curve. The differ- 
ence between the theoretical values and the oil experiment 
results can possibly be due to the fact that the viscosity 
changes in oil are large even for small temperature differences 
or due to the end effects in the measurements. 


CONCLUSIONS 


An analysis was made of the free-convection flow about a 
flat plate oriented in a direction parallel to that of the 
generating body force under the prime assumption that the 
relative temperature difference is small. 16 was found that 
the Grashof number was the principal factor determining 
the type of flow and that for large Grashof numbers the flow 
was of the boundary-layer type. The theoretical develop- 
ment was then continued to consider only the cases of large 
Grashof number because these are of most importance in 
aeronautics. 

Velocity and temperature profiles for Prandtl numbers of 
0.01, 0.72, 0.733, 1, 2, 10, 100, and 1000 were computed on 
the basis of a constant body force and plate temperature 
and agreement with experiments where the fluid was air. 
(Prandtl number of 0.72) was good. It was also demon- 
strated that velocities and Nusselt numbers of the order of 
magnitude of those obtained in forced-convection could be 
obtained in free-convection flows. 

A flow parameter and a heat-transfer parameter which are 
functions of the Prandtl number alone were derived. Calcu- 
lations of the important physical quantities such as shear 
stress, heat-transfer rate, and the like can be computed from 
these parameters. Values of the heat-transfer parameter 
obtained from an approximate theoretical development and 
from experiments compared with values computed from the 
present development showed good agreement over a wide 
range of Prandtl number (0.01 to 1000). It is shown that 
the commonly used semiempirical relation for the hent- 
transfer coefficient will yield good results only in restricted 
Prandtl number ranges. 


Lewis Еїлент PROPULSION LABORATORY 
NATIONAL ADVISORY COMMITTEE FOR ÁERONAUTICS 
CLEVELAND, Онго, October 3, 1951 


APPENDIX A 
SYMBOLS 
The following по E used in ias сЕ u |e Grashof number, 2: I. xl 
4.3, B;®, coefficients in numerical differentiation and Heo 
Qo, pm integration formulas Grx Grashof number based on X 
; g gravitational force per unit mass (or accelera- 
Cp specifie heat at constant pressure ti : 
N : . | ion due to gravity) 

F dimensionless velocity function H dimensionless temperature function 
fi components of body force per unit mass, | } heat-transfer coefficient 

i=], 2, 3 · K isothermal compressibility coefficient, 
tx negative of X-component of body force per _ E (1/p) 

unit mass P| OP | 
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k thermal-conductivity coefficient 

| characteristic length 

т, п arbitrary exponents 

N а, number, defined following equation (26) 

Nu Nusselt number, Х/К 

Ns, average Nusselt number 

P pressure 

Pr Prandtl number 

R gas constant 

r, 8, t arbitrary exponents 

T absolute temperature 

О, velocity components, 1=1, 2, 3 

Us dimensionless velocity components, 4—1, 2, 3 

и dimensionless velocity component in z-direc- 

tion с 

0 dimensionless velocity component in y-direc- 
| tion 

X, Cartesian coordinates, ?—1, 2, 3 

i dimensionless Cartesian coordinates, 4—1, 2, 3 

£ dimensionless Cartesian coordinate 

Y Cartesian coordinate 

y dimensionless Cartesian coordinate 


B coefficient of volumetric expansion, 


JE (1/p) 


ratio of specific heats 

Laplacian operator 

relative temperature difference, 8(75– Т.) 
similarity variable 

dimensionless temperature function 
step size used in numerical calculations 
absolute viscosity 

kinematic viscosity 

density 

dimensionless pressure function 

shear stress 

dimensionless density function 

stream function 


X6 09950 VER CAM PR 


$, ქ Cartesian tensor and summatjon subscripts 
8 denotes evaluation at static conditions (e=0) 


0 denotes evaluation at plate surface 


თ denotes evaluation &t undisturbed conditions 
Subscript notation is used to denote partial differentiation. 
Superscripts: 

Primes denote ordinary differentiation. 

Bars (as v or y) denote transformed dimensionless quantities. 


APPENDIX B 


NUMERICAL SOLUTION OF SIMPLIFIED BOUNDARY-VALUE PROBLEM 
By Lynn U. ALBERS 


The method is presented herein by which solutions to the 
boundary-value problem 


EP” --ЗКЕ”"—2Е?--Н=0 (B1) 


H''--3PrFH' = (B2) 


F(0) —F'(0)—0 H(0)—1 


P'(2)—H(»)—0 


were obtained for the cases of Pr equal to 0.01, 0.72, 0.733, 
1, 2, 10, 100, and 1000. "This discussion will enable the 
results to be clearly evaluated and will perhaps serve as 9 
guide in the numerical solution of similar problems. 

Each of the cases of the problem has a solution for a 
particular set of values for #”7(0) and H'(0), hereinafter 
called eigenvalues. The basic approach to the problem was 
to estimate the eigenvalues and to integrate out from zero, 
obtaining functions which satisfied equations (B1) and (B2) 
üt each step. The integration was continued until the func- 
tions 7" and H behaved in а fashion inconsistent with the 
boundary values at infinity; for example, when they became 
negative or diverged to infinity. Improved estimates of the 
eigenvalues were then made on the basis of the results of 
preceding runs and the process was repeated successively 
until a solution was obtained. 


Modifications required to overcome specific obstacles will 
be discussed after sufficient details of the basic procedure 
have been given. Then an evaluation of the accuracy of the 
numerical results will be made. 

The integration process consists of two parts, a starting 
phase and an extension phase. The starting phase begins 
with an estimate of the eigenvalues Ё””(0) and H’(0) and a 
decision on the step size x to be used. It continues with an 
iterative process of alternately computing F” and H” at 
the first four points and integrating them by five-point 
formulas. This process and that in the extension phases 
are so arranged that the differential equations are satisfied 
at each integral multiple of the step size. 

The extension phase used preceding data to integrate step 
by step beyond the fourth pomt. Diagrams of both phases 
will be given after a few preliminary explanations. 

All integration formulas used are based on the same idea. 
If a function, for example, Ё”””, 13 known at five points, there is 
а unique fourth-degree polynomial which agrees with it at 
these five points. Moreover, if the successive antideriva- 
tives (integrals) Е”, Е’, and F of F'" are known at one 
point, there are unique fifth-, sixth-, and seventh-degree 
polynomials which are successive, antiderivatives of this 
fourth-degree polynomial and which agree with #77, Е”, and 
Р, respectively, at the one point. Itis then a simple algebra 
problem to deduce from the values of F’’’ at five points and 
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Р, Е’, and F” at a single point the values of any of these 
four polynomials at any point. These results will approxi- 
mate the functions F', F", F”, and #”77 to a degree dependent 
on step size, the relative positions of the points in question, 
and the magnitude of the fifth derivative of F'"" in the 
neighborhood of these points. 

The preceding algebra problem can be и in all 
situations that arise in the starting and extension phases of 
the present problem and specific integration formulas шау 
be deduced. These formulas are discussed in the next 
paragraph. 

Let F'" be denoted at five successive points by Fo”, 
FU, FU, Fy”, and Е”. In the starting phase, these 
points are 0, x, 2x, Зк, and 4x, and Fo, В’, and PF," are also 
known. ‘Then the five sets of formulas required in the start- 
ing phase are. 


i—1, 2, 3, 4 (B3) 


4 
_ Ppr K we 
HF’ pis Au, 


. 4 
Hí-Hé р 24 447 B i=1, 2, 3, 4 (B4) 
1 = 
2 4 
Ри= Ру tisk M @ о АР" 


2 4 
Н:=Но- Но tps 2 Au BG" 


1=0 


1=1,2,3,4 (B6) 


РРР Е, + 


D. 3 в 2) As P/" 

1—1, 2, 3, 4 (B7) 
where the superscripts on the 4, and D,™ refer to the 
order of integration. 


The constants A,“ and П, may be read from the fol- 
lowing tables: | 


For Ay ара D,0: 


i=1, 9, 3, 4 (Bd) - 


For 4,49) and D,™: 


Itis now possible to diagram the steps of the starting phase 
of the integration. If each bar above a function denotes an 
improved estimate of it, and the first estimates of F,'/^, Fx”, 
Fz”, and F” are all equal to P, and similarly for tho 
H"', then the starting phase diagrams are 

(1) (В, Е’ "ru nt წესი Fy", Fe Еи), №,’ , I" 
(This diagram means that the values in parentheses are 
used with appropriate integration formulas from (B3) to 
(B7) to obtain Р, Е”, and F” at 5 —x.) 

(2) (Eh, Hy, Hy HH”, Н.”', Ну”, H,!')— Eh, Hi 

(3) (Е; FY’, XI, სს, Ни) Fy", Ну" 

(The preceding diagram means that the values in parentheses 
are substituted in the differential equations (B1) and (B2) 
to obtain Е” and H” at л=к.) 

(4) (Fo, Ру mn n. ry des E, Tt^. Е!) Е, Г, de 

(5) (Hs, Но’ P Ну’ ‚ Hi", Hs', Hs", Hy "Hh, Н; 

(6) (F5, Ey ) Е" P Da, Hx )—>Ез””", Hi” 

(7) (Е, Ee. Fy iU m T", qu Е’) 1, Ву, Py" 

(8) (Но, Но’, Но’, Ни’, Hr’, Hy", Ни") >Н, Ну 

(9) (Fs, Fy’, Fy”, Hs, Hj) Fs", Hy" 

(10) (Fo, Во’, F^ В", Р", EI, FSI el oy, FY, FAUT 

(11) (Ho, Ho’, Но”, Ну”, 2, Hs", H4) H4, Hy 

(12) (F4, XV”, FA", Н, Ви) FY", На 

It may be noted here that all four values of 7" and H” 
have been improved, and further improvement will require 
iteration of steps 1 to 12. "The start of the second iteration 


` 15 diagrammed as follows: 


(13) (Fo, d uu ur И LU 
(14) (Ho, Ну, Но’, Hi”, Н", Hs", H4") Hy Hy 

(15) (Е, Е Fy’ ) F” ) E. Hy Fy”, H” 

(16) წი В’ Fo", 72447 p, Fil", "s ––– წერა 


Successive sets of 12 steps are performed until the values of 
Еу”! and H;" no longer change. 

On. the IBM. Card-Programmed Electronie Caleulator, a 
deck of punched cards 2 inches thick sufficed to perform steps 
1 to 12. ‘Three runs of this starter deck at 3 minutes per 
run accomplished complete convergence in most cases. At 
the end of the starting process there have been computed and 
stored Fi, Fy’, Ву’, Ну, and Ну, and final estimates of 7”, 
Be. n^ Mut H”, Hy, Hy", and Ни’. 
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The extension phase has now been reached. It used a 
different set of integration formulas based on the same gen- 
eral ideas as equations (B3) to (B7). НВ”, В”, В”, 
Ез", and PF," now designate F” at any five successive 
points, and the subscript 5 denotes the next point, 


4 
PPL oes ВОВ," (88) 
. i ; 
Н = H; + 79 24 B, XI,” (89) 
к _4 
Ру РР Eng 24 Bs? Е," (B10) 
к 4 
Hy- Hor H + ey 24 Bs Н" (B11) 
/ № // I < (3) 77777 
== FO, +5 №, tow 24 B,” F (B12) 


where the Bj? and C™ are given in the following table: 


The extension phase may then be diagrammed simply аз 
follows: 

(1) (XV, XV, Jr pa dl By”, I. FE, XV”, В’ 

(2) (H,, Hi, Ну”, Hi", Ну”, Hy’, H,')H,, Hy 

(3) (Fo, Ру, Еу’, Нь, Ну) эВ", Hy" 

The values of the functions at the next point are computed 


in similar manner, where the latest sets of five values of © 


F’” and H” are used. This process advances step by step 
toward infinity. 

The extended deck of punched cards was about 3 inches 

thick and took a little over 3 minutes perrun. For Pr=0 72, 
a step size of 0.1 was used, the starting phase took 10 minutes, 
and the extension phase, about 30 minutes. When it is 
realized that about 11,000 operations were performed in the 
40 minutes per run, it may be seen that solution of the 
present problem would have been prohibitively difficult on 
desk-type calculators. Simplifications in method would have 
sacrificed accuracy or required smaller step size. 

In two-point boundary-value problems where one point is 
infinity, some problems of judgment are irivolved as to where 
infinity is, and as to when a satisfactory approximation to a 
solution has been obtained. Ih most cases this question was 
settled for the present problem by calling a run satisfactory 
when it fell between two runs for which Æ” and H did not 
differ at important points in the fourth decimal place, and 


for which Æ” and Н flattened out at zero, correct to four 
decimal places. 

Certain difficulties were met in the attempt to use the 
basic procedure previously discussed. "These necessitated 
certain modifications. 

For Pr=2, 10, 100, and 1000, H would settle down to 
zero at an early stage; but while Ё” was still coming down, 
H” would begin to oscillate and these oscillations increased 
and fed back into all other functions. This trouble was 
avoided by the following modifications: It is a consequence 
of equation (B2) that 


H’ (n)=H' (0) exp (—3Pr Í Pdi) 
нер (—3Pr | rou (B13) 


The extension phase was modified to require the additional 


- Integration formulas 


შე (814) 


720 £i 


5 
HH os 2-1 А.Н; (В15) 
where 4,— —19, 4:=106, A3— — 264, 4,—646, and A4—251. 

These formulas were used along with equations (B8) to 
(B10) according to the following diagram: 

(1) (Fa Еи, FU, EY”, В", В", Ру", Bl!" Fs, E, FU 

` (2) Gh, Fs, #3, Fa, Fs, Hi’) -H’ by means of (B14) and 
(B13) 

(3) (Hr, Ну, Ну, Ни, Ну, Н)—>Нь by means of (B15) 

(4) (Fs, Ру, Fo", Hy) Be” 

The value #%””” is discarded апа F’”’ at the last five points 
is used to repeat the whole process again and again ad 
infinitum. As long as # stays positive, H’ is guaranteed to 
approach zero and A will flatten out to some value and not 
oscillate. 

For Pr=0.01, 0.72, 0.733, and 1, Ње Е””’ began to oscillate 
at an advanced point and these oscillations grew and fed into 
the other functions. For all cases but Pr=0.01, the oscilla- 


tions appeared very late, near the end of the run, and a 


suitable halving of step size when oscillation was detected 
in the fourth differences of F’’’ was sufficient to avoid the 
difficulty. But for the 0.01 case, oscillations of F’’’ appeared 
early in the run, namely, soon after the peak in F. These 
oscillations were found to be step-size connected, so that 
reduction of the step to 0.02 avoided them. Even then 
oscillations in Ё”” would begin to appear every 25 steps or 
so, and these were smoothed out regularly by repeated runs 
of a deck similar to the starter deck. Each run under these 
conditions took about 16 hours, ЗНОВ this the most diffi- 
cult case to solve. 
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APPENDIX C 
DERIVATION OF FLOW PARAMETER 


By definition the shear stress is given by 


= (=>) 
T= H0 OY о 


To express (QU/O Y), in terms of the known function XVI), 
use can be made of equations (60) and (62). Then 


(C1) 


Substitution of this expression into equation (C1) yields the 
flow parameter 
ONE UIS =F” (0) 
(&Grx9)' (va uo X?) 


Note that from the general derivation, the flow parameter 
contains the viscosity evaluated at two different points. 
Recall, however, that the analysis has shown that to a first 
approximation the variation of viscosity with temperature 
can be neglected. Thus the viscosity can be taken as 
constant in the entire flow field. 


APPENDIX D 
DERIVATION OF HEAT-TRANSFER PARAMETER 


The local Nusselt number is defined as 


„გ ...-.-. 


То express (0 7/0Y), in terms of the known function H(y), 
use is made of equations (61) and (62). Thus 


afta (Em 


Substitution of this expression into equation (D1) yields 
the heat-transfer parameter 


‘of this parameter. 


The heat-transfer parameter as given by equation (D2) 
is, аз was previously stated, a local parameter. It is often 
desired to compute the average (over the length X) value 
To this end, the Nusselt number (as 
given in equation (D1)) must be defined in terms of an 
average heat-transfer coefficient and the quantity thus ob- 
tained must then be integrated over the length X and 
divided by X. This procedure yields the result 


=> (VI 


Nu =H’ (0) (D2) | It is from this last equation that the 0.75 reduction factor 
(Grx/4) ! previously discussed was obtained. 
REFERENCES 


1. Eckert, E. R. G., and Soehngen, E. E.: Studies on Heat Transfer 
in Laminar Free Convection with the Zehnder-Mach Interferom- 
eter. Tech. Вер. No. 5747, A. T. I. No. 44580, Air Materiel 
Command (Dayton, Ohio), Dec. 27, 1948. 

2. Schmidt, Ernst, und Beckmann, Wilhelm: Das Temperatur- und 
Geschwindigkeitsfeld vor einer Würme abgebenden senkrechter 
Platte bei natürlicher Konvektion. Tech. Mech. и. Thermo- 
dynamik, Bd. 1, Nr. 10, Okt. 1930, pp. 341-349; cont., Bd. 1, 
Nr. 11, Nov. 1930, pp. 391-406. 

3. Schuh, H.: Boundary Layers of Temperature. Reps. & Trans. 
1007, AVA Monographs, British M. А. P., April 15, 1948. - 

4. Ostrach, Simon: А Boundary Layer Problem in the Theory of Free 

; Convection. Doctoral Thesis, Brown Univ., Aug. 1950. 

5. Lewis, J. А.: Boundary Layer in Compressible Fluid. Мопо- 
graph V, Tech. Вер. No. F-TR-1179-ND (GDAM A-9-M V), 
Air Materiel Command (Dayton, Ohio), Feb. 1948. (Analysis 
Div., Intelligence Dept. Contract W33-038-ac-15004 (16351) 
with Brown Univ.) 

0. Zemansky, Mark W.: Heat and Thermodynamics. Second ed., 
McGraw-Hill Book Co., Inc., 1943, pp. 26-27. 


7. Eckert, E. В. G.: Introduction to the Transfer of Heat and Mass. 
MeGraw-Hill Book Co., Inc., 1950. 
8. Lagerstrom, Paco, Cole, Julian D., and Trilling, Leon: Problems in 
the Theory of Viscous Compressible Fluids. Guggenheim Aero. 
Lab., C. I. T., March 1949, pp. 38-40. (Office of Naval Res. 
Contract N6onr—Task VIII.) í 
9. Friedrichs, К. O., and Wasow, W. R.: Singular Perturbations of 
Non-Linear Oscillations. Duke Math. Jour., vol. 13, 1946, pp. 
367—381. 
10. Eckert, Е. Б. G., and Jackson, Thomas W.: Analysis of Turbulent 
Free-Convection Boundary Layer on Flat Plate. NACA Rop. 
1015, 1951. (Бирегведез МАСА ТМ 2207.) 
11. McAdams, William H.: Heat Transmission. Second ed., MoGraw- 
Hill Book Co., Inc., 1942, p. 242. 
12. Lorenz, Hans H.: Die Würmeübertragung von einer ebenen, 
senkrechten Platte an: Ol bei natürlicher Konvektion. Zs. f. 
tech. Phys., Nr. 9, Heft 15, 1934, pp. 302-360. 
18. Saunders, О. A.: Natural convection in liquids. Ргос. Roy Sos. 
(London), ser. À, vol. 172, no. 948, July 19, 1939, pp. 55-71. 


FREE-CONVECTION FLOW AND HEAT TRANSFER ABOUT А FLAT PLATE 


TABLE I—FUNCTIONS F AND Н AND DERIVATIVES FOR VARIOUS PRANDTL NUMBERS 
(8) Prandti number, 0.01 


Е" Н Ir Е F” H H’ 

0. 9882 | 1.0000 | —0. 0812 4.16 I. 9319 0.4037 | —0.0615 | 0.6741 | —0.0723 
ol ‚ 0048 8 . 8868 . 9919 —. 0812 4.18 1. 9400 . 4024 —. 0614 6726 —. 0722 
.2 .0184 .1774 . 7891 . 8838 —. 0812 4.20 1. 9180 .4012 —. . 6712 —. 0721 
.9 0309 . 2616 . 6942 ‚ 9756 —, 0812 4.24 1.9640 . 9088 — „6833 —. 0720 
.4 0684 . 9104 . 6030 . 9675 —.0811 4.23 1. 9709 . 3004 —. 0005 „ 6654 —. 0718 
. 5 .1029 » . 5163 . 0504 —. 0811 4. 32 1. 0957 . 3839 —. . 6620 —.0716 
.6 1426 „ 4108 . 4349 „9513 —. 0811 4. 36 2. 0114 . 9015 — . 6597 —. 0714 
7 . 1866 > „ 3595 . 0432 —.0811 4. 40 2. 0270 . 3891 —. 0595 . 6568 —. 0713 
.8 » 2242 „4919 . 2006 . 9351 —. 0810 4. 50 2. 0657 . 9832 —. 0587 . 6497 —. 0708 
-9 . 2848 . 5178 . 2235 . 0270 —. 0809 4. 60 2. 1037 8774 —. 0579 6427 —. 0704 
1.0 - 3376 . 5379 . 1784 „4189 —. 0809 4, 70 2.1411 3716 —. 0572 . 6357 —. 0699 
1.1 . 3922 . 5527 . 1253 . 9108 —. 0808 4. 80 2. 1780 3680 —. 0564 . 6287 —. 0695 
1.2 . 4480 . 5831 . 0839 . 9028 —. 0807 4.90 2. 2143 3004 —,. 0556 6218 —. 0690 
1,3 . 5047 . 5697 . 0489 . 8047 —. 0800 5.00 2. 2501 3548 —. 0549 6149 —. 0696 
1.4 . 5818 . 5781 . 0198 . 8866 —. 0804 5.10 2. 2853 3494 —, 0641 . 6081 —. 0631 
1.5 ‚ 0102 - 5739 — . 8788 —. 0803 5. 20 2. 3200 3440 —. 0534 6013 —. 0676 
1.6 . 6765 . 5725 —, 0229 . 8706 —.0801 5. 40 2. 8877 3335 —. 0520 „5578 —. 0667 
1.7 . 7338 . 5604 —. 0379 . 8620 —. 0800 5.60 2. 4534 3232 —. 0508 . 5746 —. 0657 
1.8 . 7904 . 5850 —. 0193 . 8546 —. 0798 5.80 2. 5170 . 3132 —. 0493 5815 —. 0818 
1.9 . 8466 . 5596 —. 0579 . 8466 —. (7 6. 00 2. 5787 $035 —. 0479 . 5487 —. 0638 
2.0 . 9023 . 5535 —. 0041 . 8387 —. 0794 6. 20 2. 6384 . 2040 — . 5360 —. 0628 
2.1 ‚ 9573 . 5469 —. 0685 . 8307 —. 0792 6. 40 2. 6963 . 2840 —. 0454 . 5236 —. 0518 
2.2 1. 0117 „ 5399 —. 0714 .8228 —. 0789 6. 60 2. 7524 . 2759 — . 5113 —. 0608 
2.3 1. . 5320 —. 0; . 8150 —. 0787 8. 80 2. 8067 . 2671 —. 0430 „4993 —. 0593 
2.4 1.1182 . 5253 — 07 . 8071 —. 0784 7. 00 2. 8593 . 2586 —. 0419 . 4874 —. 0588 
2.5 1.1703 . 5178 —. 0748 „7993 —. 0782 7.40 2. 9594 2423 — 0397 „4643 —. 0568 
2.6 1.2217 . 5104 —. 0745 7915 —. 0779 7.80 3. 0532 . 2260 — „ 4420 —. 0517 
2.7 1. 2724 „ 5029 —. 0741 7837 —. 0776 8. 20 3. 1411 . 2123 —. 0356 . 4205 —. 0527 
2.8 3223 „4956 —. 0735 7760 —. 0773 8. 60 3. 2232 . 1984 —. 0337 . 3998 — 0508 
2.9 1. 3715 „4852 —.0728 . 7052 —. 0770 9. 00 3. 2099 1853 —. 0319 .3799 —. 0489 
3. 0 1. 4200 .4810 —.0719 . 7006 —. 0766 9. 40 3. 3715 .1729 — „3607 —. 0470 
3.1 1. 4077 .4739 —. 0711 . 7529 —. 0763 9. 80 3. 4383 . 1612 —. 0285 . 3423 —. 04151 
3.2 1. 5148 $ —. 0701 . 7453 —. 0760 10.20 3. 5005 1501 —. 0270 .3246 — 
3.3 1. 5611 „4595 — .7377 —. 0756 10. 60 3. 5584 . 1396 —. 0254 . 3069 —. 0413 
3.4 1. 6067 . 4530 —. 0683 . 7302 —. 0753 11.00 3. 6123 . 1297 —. 0241 . 2915 —. 0397 
3. 45 1. 6293 . 4406 —.0078 . 7204 —. 0761 11, 40 3. 6822 1203 = . 2759 —. 0380 
3.5 1. 6517 . 4462 —. 0874 . 7227 —. 0749 11.80 3. 7086 .1114 —. (215 . 2810 —. 0364 
ვ. 55 1. 6739 . 4428 — „7189 —. 0747 12. 20 3. 7514 1031 — . 2468 —. 0348 
3.6 1. 6960 „ 4395 =, . 7162 —. 0745 12. 60 3. 7911 . 0952 —. 0192 . 2332 —. 0332 
3. 65 1.7179 .4362 —, 0660 .7115 —, 0743 13. 00 3. 8276 0377 —. 0181 . 2202 —.0317 
3.7 1. 7396 „4329 —, . 7078 —, 0741 13. 60 3.8771 . 0773 —.0165 . 2018 —. 0298 
8. 75 1. 7611 . 4296 —. 0851 . 7041 —. 0739 14. 20 3. 9206 0070 —. 0151 . 1847 —. 0276 
3.8 1. 7825 „ 4201 —. . 7004 —. 0738 14.80 3. 9587 0592 —. 0138 .1 —. 0257 
4. 85 1. 8308 .4232 —. 0642 . 6067 —.0738 15. 40 3. 9918 . 0513 —. 0126 ok —. 0239 
3.9 1.8249 . 4200 —. 0038 . 6030 —. 0734 16. 00 4. 0201 . 0441 —. 0114 xL —. 0223 
3. 95 1.8453 . 4168 —, 0633 . 6893 —. 0732 17. 00 4. 0591 . 0335 —. 0097 1189 —.0197 
4.00 1.8685 . 4136 —. 0829 . 0857 —. 0729 18. 00 4. 0880 . 0246 —. 1003 —. 0175 
4.05 1. 8871 . 4105 —. 0025 . 0821 —. 0727 19. 00 4.1088 0171 —. 0069 . 0839 —. 0154 
4.10 1. 9076 .4074 —. 0520 . 0784 —. 0725 20. 00 4. 1226 0108 = . 0604 —.0137 
4.12 1. 9157 . 4061 —. 0619 6770 —. 0725 21.00 4. 1308 3 —. 0046 . 0565 —. 0121 
4.14 1, 9238 „ 4049 —. 0617 6755 —. 0724 22.00 4. 1343 0015 = . 0452 —. 0107 

(b) Prandtl number, 0.72 (с) Prandtl number, 0.733 


ol el —. 5079 
.2 .2 —. 5071 
3 .3 —. 5050 
.4 ‚4 —. 5012 
5 5 —. 4953 
0 8 —. 4870 
-T -T —. 4783 
„8 .8 —.4632 
.9 .9 —.4478 

1.0 1.0 = 
1.1 1.1 —.4110 
1.2 1.2 —.3902 
1.3 1.3 —. 3684 
Li. 1.4 —: 3458 
1.5 . 1.5 —.3228 
1.8 . 3490 1.6 —. 2998 
1.7 . 9707 1.7 —. 2771 
1.8 .3910 1.8 —. 2549 
1.9 . 4100 1.9 —.2335 
2.0 „4277 2.0 —. 2130 
2.1 „4440 2.1 –.1937 
2.2 „ 4690 2.2 —.1754 
23 .4728 2.3 —. 158i 
2.4 4854 2.4 —.1427 
2, 5 .4969 2.5 —. 1281 
2.6 . 5074 2.6 —. —. 1148 
2.7 . 5168 2.7 —. . 0843 | —. 1026 
2.8 . 5254 2.8 -—. .0746 | —.0916 
2.9 . 5332 2.9 =, .0650 | —. 0816 

3.0 . 5401 3.0 —. .0583 | —. 
3.1 . 5464 8. 1 —. .0514 | —.08014 
3.2 .5520 3.2 -—. 0454 | —.0571 
3.3 . 5571 3.3 а 0400 | —.0506 
3,4 . 5816 3.4 .0352 | —.0448 
3. 6 . 5692 3. 6 0273 | —.0350 
3.8 . 5735 3.8 0211 | —.0272 
4.0 5801 4.0 .0163 | —.0211 
4.2 „ 5338 4.2 „5307 .0126 | –.0161 
4.4 . 5888 4.4 . 5837 .0097 | —.0127 
4.6 . 5391 3.6 . 5860 .0074 | —.0098 
4.8 . 5908 4.8 ‚5879 .0057 | —.0076 
5.0 „ 5922 5.0 . 5893 .0044 | —. 005$ 
6.2 .5933 5.2 . 5905 .0033 | —. 0045 
5.4 . 5941 5.4 .5914 .0025 | —. 0035 
5.8 .5951 5.6 . 5921 .0019 | —. 0027 
6. 2 .5957 5.8 . 5927 .0015 | —.0021 
6.8 . 5960 6.0 . 5932 „0011 | —.0016 

7.3 . 5961 6. 4 . 5938 .0008 | —. 

6.8 . 5913 . 0003 —. 

7.2 .5946 . 0001 | —. 

7.6 . 5949 . 0000 —. 
8.0 . 5951 .0000 | —. 0001 
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TABLE L—FUNCTIONS F AND H AND DERIVATIVES FOR VARIOUS PRANDTL NUMBERS—Continued 


Е 
(d) Prandtl number, 1 (e) Prandti number, 2 
| - 
2 F F Fv H H' 
0 0. 0000 0. 0000 0. 5713 1.0000 | —0. 7165 
^ -1 2 „0523 „ 4749 „9284 —. 7161 
.2 0101 „0952 . 3881 8569 –, 
.8 . 0215 .1297 . 3049 7858 —. 7069 
М .4 0358 «1585 . 2318 7157 —. 6919 
а .5 0525 „1763 . 1668 6470 —. 6763 
а .6 0709 . 1901 . 1095 é —. 6523 
გ -7 „0904 „1985 . 0602 „6168 —. 6215 
გ 8 . 1104 . 2024 . 0185 .4 — 
.9 . 1307 . 2024 | —.0161 —. 5443 
1.0 ~ . 1508 .1991 | —.0440 . 3476 —. 5002 
1.1 „1705 .1938 | —*0659 к —.4543 
1.2 2 „1884 —. 0824 — 
1.3 2077 .1775 | —.0942 2183 —. 3619 
1.4 . 2250 .1677 | —.1020 1844 —. 31 
1.5 . 2412 .1572 —. 1063 „1547 —. 2763 
1.6 . 2564 . 1465 —. 1078 . = 
1.7 . 2706 .1357 | —.1071 . 1069 —.2032 
1.8 . 2836 .1252 | —.1046 . 0882 —. 1721 
1.9 . 2956 .1149 | —.1008 . 0724 —. 1446 
2.0 „3066 21 —. 0960 . 0592 —.1207 
2.1 . 3166 .0957 | —.0906 . 0482 —. 1001 
2.2 . 3257 „0589 | —.0849 . 0391 — 
2.3 . 3340 .0787 | —.0789 . 0316 —. 0677 
2.4 . 9415 . 0711 —. 0729 . 0254 —.0553 
2.5 . 9483 .0041 | —. 0671 . 0204 —. 0450 
2.6 . 9543 . 0577 —.0614 . 0164 —. 0364 
2.1 2.7 . 9508 .0518 | —.0500 . 0131 —. 0204 
2.8 2.8 . 3647 " —. 0509 „0105 —. 0237 
2.9 2.9 . 3001 .0417 | —.0461 . 0083 —. 0190 
3.0 3.0 „3731 .0373 | —.0416 . 0066 —. 0152 
3.1 3.1 . 9760 . 0333 | —. 0375 0053 —, 0121 
3.2 3.2 . 3798 .Q208 | —.0338 0042 — 
3.3 3.3 . 9826 .0268 | —. 0303 0033 —. 0077 
34 3.4 . 3851 —.0272 0026 —.0061 
3.6 3.6 „2893 . 0188 | —.0218 0017 —. 0038 
3.8 3.8 3927 .0140 | —.0174 0010 —. 0024 
4.0 4.0 . 9933 .0018 | —. 0138 . 0007 —.0015 
4.2 4.2 3974 . 0091 —. 0109 . 0004 — 
4.4 4.4 3991 .0074 | –.0088 . 0003 —. 0006 
4.6 4.6 . 0059 —. 0068 . 0002 — 
4.8 4.8 . 4015 .0047 | —.0054 . 0001 —. 
5.0 5.0 .4023 Py —. 0042 . 0001 — 0001 
5.2 5.4 „4035 * — . 0001 —, 0001 
5.5 5.8 .4043 .0015 | —.0016 . 0000 . 0000 
6.0 6.4 „ 4049 2 —. 0008 . 0000 . 0000 
6.25 7.0 .4053 . 0005 | —.0001 . 0000 . 0000 
8.0 . 4036 .0002 | —.0001 . 0000 . 0000 
9.0 „4058 „0001 „0000 „0000 . 0000 
10.0 . 4050 „0001 . 0000 . 0000 . 0000 
11.0 . 4059 . 0001 . 0000 . 0000 . 0000 


3 F Е Е" H H’ 9 F F’ F” H H’ 

0 0. 0000 0. 0000 0.4192 1.0000 | —1.16 . 0. 2517 1.0000 | —2.191 
.I . 0019 . 0371 .3251 .8832 | —L 1671 . . 0001 . 0000 . 2274 .0462 | —2.191 
.2 2 . 0654 . 2428 .4070 | —1.1521 . 050 . 0003 .0114 . 2044 .8905 | —2.188 
.3 . 0147 . 0881 . 1723 .6534 | —1.1155 |" - 075 . 0006 . 0162 . 1828 .8359 | —2. 180 
.4 . 0241 . 1003 . 1134 6148 | —1.0526 . 100 . 0011 Р .1626 «7815 | — 2.166 
.5 . 0346 . 1091 . 0655 „4437 —. . 125 0017 . 0244 . 1438 .7276 | —2.144 
.6 . 0458 . 1137 . 0279 .3527 —. 8515 .150 - 0023 . 0277 . 1203 . 6744 | —2.113 
-T " .1150 .| —.0011 . 2133 =, . 175 . 0031 . 0307 . 1101 . 6221 | —2.071 
.8 .0687 . —.0221 „ 2004 —. . 200 . 0039 . 0832 . 5700 —2. 018 
.9 А .1107 | —. 0367 . 1519 —. 4849 . 225 . 0047 . 0355 . 0816 .5213 | —1.954 
1.0 . 0908 -1066 | —.0462 . 1090 —.3753 . 250 . 0056 . 0373 . 0692 .4733 | —1.880 
1.1 „1012 .1016 | —. . 0763 —. . 875 . 0066 . 0389 . 0580 .4273 | —1.786 
1.2 «1111 $ —. 0545 .0522` | —.2045 .300 . 0076 . 0402 . 0479 .3836 | —1.704 
1.3 " .0908 | —. . 0349 —. 1445 .325 . 0086 . 0413 . 0389 .3422 | —1. 604 
1.4 1293 .0853 | —. . 0228 —. 0993 350 . 0006 . 0422 . 0309 .3034 | —1.498 
1.5 . 1376 .0799 | —. . 0146 —. 0665 |- .375 .0107 .0429 0238 .2674 | —1.383 
1.6 - 1453 .0748 | —. i —. 0435 -400 0118 . 0434 0176 .2341 | —L 276 
1.7 . 1525 .0699 | —. . 0056 —. .425 0129 . 0438 0123 .2036 | —1.163 
1.8 . 1593 .0852 | —. 0453 . 0034 —. 0174 „460 0140 „0140 0076 .1759 | —1.052 
1.9 . 1656 . 0508 | —. .0020 —. 0107 „475 „0151 . 0442 . 0036 . 1510 —.9430 
2.0 . 1714 .0587 | —. . 0012 —. 0065 . 500 . 0162 . 0442 . 0002 . 1287 —. 8383 
2.1 . 1769 .0538 | —.0375 " —.0038 . 525 . 0173 .0442 | —.0026 . 1089 —. 7404 
2.2 . 1820 .0491 | —. 0351 . 0004 —. 0022 . 550 . 0184 .0441 | —.0050 . 0916 —. 0178 
2.3 . 1888 .0458 | —. . 0002 —. 0013 . 575 . 0195 . 0439 | —.0070 . 0765 —. 5021 

2.4 . 1912 . 0426 | —. - 0001 —.0007 . 600 . 0206 .0437 | — . 0634 —. 4837 

2.5 . 0396 | —. . 0001 —. 0004 . 65 . 0227 .0432 | –.0111 . 0427 —. 3496 
2.6 . 1991 .0369 | —. . 0000 —. 0002 . 70 . 0219 .0426 | —. 0128 . 0280 —. 240 
2.7 - 2027 . 0343 | —. . 0000 —. 0001 .75. . 0270 .0420 | —. 0135 . 0178 —. 1657 
2.8 . 2060 0319 | —. . 0000 —. 0001 . 80 . 0291 .0413 | —.0140 0111 —. 1088 

2.9 . 2090 .0297 | —.0216 . 0000 . 0000 .85 . 0311 . 0406 | —. 0142 . 0067 —. 0603 

3.0 . 2119 .0276 | —.0201 . 0000 . 0000 . 90 .0331 .0399 | —.0142 s —. 0428 

3.1 . 2146 . 0256 | —. 0187 . 0000 . 0000 . 95 . 0351 .0392 | —.0141 . 0022 — 

3.2 . 2110 .0238 | —. 0174 "e „0000 1.00 . 0371 .0385 | —.0140 . 0012 —. 0149 

3.3 . .0221 | —.0162. $ - 0000 L 10 . 0408 . 0371 | —.0136 . 0004 = 

3.4 . 2215 .0208 | —. 0151 . 0000 . 0000 1.20 . 0445 .0357 | —. 0132 . 0001 —. 0013 

3.6 . .0178 | —.0131 . 0000 . 0000 1.30 . 0480 .0344 | —.0128 . 0000 —. 0003 

3.8 . 2286 .0153 | —. 0113 . 0000 . 0000 1.40 . 0514 .0332 | —. 01% . 0000 —. 0001 

4.0 .2314 .0132 | —. . 0000 . 0000 L 50 . 0546 .0320 | —.0119 . 0000 0000 

4.2 » .0114 | —. . 0000 0000 1. 60 . 0578 .0308 | —.0116 . 0000 0000 
4.4 2360 .0098 | —. . 0000 . 0000 1, 70 . 0608 .0297 | —.0112 . 0000 0000 
4.6 .2319 .0085 | —. . 0000 . 0000 1.80 . 0637 .0286 | —.0108 . 0000 0000 

4.8 . 2891 . 0073 | —. . 0000 . 0000 1.90 . 0065 .0275 | —.0104 . 0000 0000 
5.0 . 2408 .0063 | —.0047 . 0000 . 0000 2.0 . 0692 . 0285 | —.0101 . 0000 0000 
5.4 . 2430 .0047 | —. . 0000 . 0000 2.1 . 0718 .0255 | —.0097 . 0000 0000 
5.8 . 2446 .0035 | —. . 0000 . 0000 22 . 0743 .0245 | —.0094 . 0000 0000 
6.2 . 2458 .0026 | —.0019 . 0000 . 0000 2.3 . 0767 .0236 | —.0091 0000 0000 
7.0 . 2174 .0014 | —. 0011 . 0000 . 0000 2.4 . 0790 .0227 | —. 0088 . 0000 4 
8.0 n .0007 | —. 0005 . 0000 . 0000 2.6 Я .0210 | —. 0082 0000 0000 
9.0 . 2489 .0003 | —.0002 . 0000 . 0000 2.8 .0874 .0195 | —.0076 0000 0000 

10.0 . 291 .0002 | — 0001 . 0000 . 0000 3.0 . 0912 .0180 | —.0071 . 0000 . 0000 
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TABLE I—FUNCTIONS F AND Н AND DERIVATIVES FOR VARIOUS PRANDTL NUMBERS—Concluded 
(g) Prandtl] number, 100—Concluded (h) Prandtl number, 1000 
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